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TECHNICAL  DISCUSSION 


The  grant  period  was  devoted  to  a study  of  two  classes  of  non- 
linear optical  processes: 

(1)  Nonlinear  optical  mixing  of  microwave  radiation  at  4 GHz, 
with  the  944  cm  ^ P(20)  line  of  CO^  in  a NH2D  gas.  This  analysis  formed 
the  basis  for  a series  of  experiments  at  the  Hughes  Research  Laboratories, 
Malibu,  where  the  effect  was  observed  and  agreement  with  the  theory  es- 
tablished. 

The  work  is  described  in  the  enclosed  reprint  (Appendix  1). 

(2)  Conjugate  four-wave  mixing  for  time  reversed  propagation:  A 
theoretical  analysis  of  four-wave  mixing  in  oases  was  conducted.  Special 
emphasis  is  on  the  use  of  this  effect  for  correcting  for  aberrations  in 
optical  propagation. 
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APPENDIX  1 


Stark-Induced  Three-Wave  Mixing  in  Molecular 

Gases — Part  I:  Theory 

RICHARD  L.  ABRAMS,  mumblr  ieee,  AMNON  YARIV,  fellow,  ilee,  and  PO  CHI  A.  YEH 


/I6ifrar/- Application  of  a dc  electric  field  to  a gaseous  system 
destroys  the  basic  inversion  symmetry  and  allows  lliree-wave  miving 
processes  to  occur.  A Lbcorclical  derivation  of  tins  effect  under  con- 
ditions of  resonantly  enhanced  nonlincaritics  is  given  for  a three-level 
lystcni.  Calculations  are  presented  for  mixing  of  a CO;  laser  with 
X-Gllz  iiiictowaves  in  the  molecule  Nil;!),  producing  single  lower 
sideband  ladiation. 


Manuscript  received  September  2,  1976  This  work  was  supported  in 
part  by  the  Advanced  Research  I'rojects  Agency,  monitored  by  the 
Office  of  Naval  Ucseatcli,  and  in  part  by  the  Army  Research  Office, 
Duriiam,  NC. 

K I Abrams  is  with  Hughes  Kesearcli  Laboratories,  Malibu,  CA 
9026.V 

A.  Yariv  and  P Ych  are  with  the  California  Institute  of  Technology, 
Pasadena.  CA  91 109. 


I.  Introduction 

'J^T ONLINEAR  optical  inixing  in  atomic  vapors  has  been 
i.  xj  demonstrated  for  a number  of  different  processes  in- 
cluding third  harmonic  generation  (l]-[4],  dc-induced 
second-harmonic  generation  (5|,  (6],  infrared  upconversion 
[7],  (8),  and  mullii'hoton  generation  of  new  wavelengths 
(9) , Resonant  enhancement  and  phase  matching  of  three- 
photon  processes  has  led  to  rather  impressive  conversion  ef- 
ficiencies for  certain  inietactions  in  atomic  vapors  [2) -[4]. 
In  this  paper  wc  discuss  Ihiec-wave  mixing  processes  in  mole- 
cules whcic  resonant  enhancement  is  aclitevcd  via  Staik  tuning 
of  the  moicculai  cneigy  levels.  This  interaction  and  its  sub- 
sequent experimental  observation  (10)  suggest  a new  type  of 
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clcctiooptical  eflcct,  namely  sinple  siilcband  pcneration  by 
applied  mietowave  Ireqiiencies  A theoieiical  derivaiiDti  and 
calculations  of  the  interaction  for  a three-level  system  are 
presented  liere,  specialized  to  (he  case  ot  a particular  molecule 
(N'lljD).  The  experitnental  observations  are  discussed  in  the 
following  paper  (10) . 


11.  Theoky 


The  application  of  a dc  electric  field  to  a gas  introduces  a 
preferred  spatial  direction,  thus  destroying  the  inversion  sym- 
metry. The  second-order  induced  polariz.ation  amplitude  can 
then  be  related  to  the  product  of  tlie  field  (complc.v)  ampli- 
tudes by 


' Q 


= J 


(j,  =ca,-ta 
aOy 


(1) 


Choosing  the  direction  of  the  dc  field  as  2,  the  allowed  JaPr 
are  cfjj. , o'.,-,,  and  cf,-,,,  where  i = .r  or  y. 

In  searching  for  a candidate  gas  in  which  to  observe  the 
effect,  one  should  look  for:  1)  molecules  with  a strong  perma- 
nent dipole  moment  or  2)  molecules  which,  in  the  presence  of 
a dc  field,  acquire  a large  dipole  moment  so  that  the  presence 
of  the  dc  field  constitutes  an  appreciable  perturbation. 

A molecule' meeting  criterion  2)  is  NHjU.  The  molecule  has, 
among  others,  the  three  levels  shown  in  Fig.  1,  which  can  be 
Stark-tuned  into  simultaneous  resonance  with  the  /^(20)  line 
of  the  CO;  laser  [lll-[14j  and  microwave  radiation  near  4 
Cllz.  as  shown.  Inis  should  lead  to  a strong  resonant  niixing 
of  the  f’(20)  line  (of  frequency  <03 /2n)  and  the  microwave 
field  at  coj/27r  = 4 GHz,  giving  rise  to  the  difference  frequency 
radiation  at  oz,  =co3  - co,  when  the  Stark  field  is  near  - 
3570  V/cni.  Levels  1 and  2 belong  to  the  lowest  vibrational 
stale  (I'j  = 0)  and  have  molecular  angular  momentum  quantum 
numbers  7 = 4 and  1,A/|  = 4.  The  subscripts  04  and  14  corre- 
spond to  the  standard  asymmetric  top  designation  [15],  The 
symbols  a (asymmetric)  and  s (syminetiic)  refer  to  the  parity 
of  the  inversion-split  sabrational  wave  functions.  The  applica- 
tion of  an  electric  field  Ep^  causes  an  admixture  of  the  wave 
functions  l-toau)  and  |4,4S>  which  is  due  to  a noiivanishing 
matrix  element  of  the  molecular  dipole  operator  connecting 
the  two  states.  This  ailmLxture.  whicii  will  soon  be  shown  to 
be  responsible  for  the  nonlinear  mixing,  disappears  at  zero  dc 
field.  The  paiameier  A appearing  in  the  expression  for  the 
wave  functions  corresponds  to  the  energy  splitting  /f;  ' 
between  the  two  low-lying  states  and  is  given  by 

A=  lM<‘io^alyJ.4,^s>l^(Ep,y  (2) 

while  the  admixture  w'ave  functions  are 


11)  = ^ |x/r+67A14o4a)  + VT^6/A|4,4S>], 


|2> 


= ^j[\/l  |4oaa>-  |/l  +^i4,4S> 


where  6 is  the  zc  -tting  and  p_.  is  the  projection  of 

the  molecular  dip  ■ operator  along  the  direction  of 

the  dc  field. 

The  expression  for  the  nonlinear  dipole  moment  of  an 
NlljO  molecule  depends  on  matrix  elements  which  can  be 
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I ig.  1 . Sonic  of  the  energy  levels  lelevant  to  the  derivation  of  the  non- 
linear coefficient. 


determined  from  linear  absorption  data  as  well  as  from  the 
data  on  Stark  splitting.  This  makes  possible,  in  principle,  a 
precise  theoretical  derivation  of  the  nonlinear  mixing  behavior 
of  this  molecule  and  of  its  parametric  dependencies. 

Applying  second-order  perturbation  theory  (16,  p.  556]  to 
the  three-level  system  of  Fig.  l,and  keeping  only  the  resonant 
(i.e.,  with  near  vanishing  denominator)  term,  leads  to  the  fol- 
lowing expression  tor  ihe  poiarization  generated  at  co,  = 003  - 
coj  by  the  applied  fields  at  403,003; 

D(2)f  \ _ ^ I ^1(7*  '£^3)ia(M'f:?)3i  (^0)3* 

“ 47?^  1 (I'n  + *(‘*^3  ■ W31 )]  [r32  + i(tO|  - 0)33)]  . 

_ (A';  - 1 

(r  12  + 1(403  ■ <^21  )l  (J  32  f(Wi  - 4*332)1  J 

• exp  (tcoi  r) -h  c.c.  (3) 

where  A',-  is  the  population  density  of  level  1 with  Ej  - Ej  = 0. 
At  thermal  equilibrium  A'j  — A'j,  and  the  main  contribution 
toFi'>  is  from  the  first  term,  the  one  proportional  to  A'j . 

At  zero  dc  field  the  matrix  clement  (Ma)32  is  zero.  Tltis  is 
due  to  the  fact  that,  as  can  be  shown  by  group  tlieoictic 
arguments,  only  the  molecular  dipole  moment  along  tlic  b of 
NII2D  axis  (71ft)  may  possess  a nonvanishing  matrix  clement 
<5o5fl|Pftl4i4S>,  but  Mft  =0  due  to  the  basal  plane  symmetry 
of  NHjD.  It  follows  from  (3)  that  for  E^^  =0  no  frequency 
mixing  takes  place.  Wlien  E^^  ^O  the  ground  state  wave 
function  |4o4iz>  is  admixed  into  level  2 as  shown  in  Fig.  I. 
This  results  in  a nonvanishing  matiix  element  (714)32  pro- 
portional to  (Sosal/a^  |4|J4J>. 

I'or  £'2112,  a = x,  and  t'sll-t  we  find,  using  the  admixed 
wave  functions,  that  the  triple  matrix  clement  product  ap- 
pearing in  (3)  is  given  by 

= -471c  <4o4i7l7ixl5osi»>* 

•£dc4f*(fr)-  W 

The  dependence  of  the  triple  matrix  clement  product  on  the 
dc  clectiic  field  IS  contained  in  the  factor  Ej^  f>IS^  with  A Ihe 
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energy  separation  between  levels  2 and  1 , as  given  by  (2).  The 
nonlinear  mixing  is  thus  absent,  t.e.,  = at  zero  field 

{t'ac  ~ 0)  and  at  very  high  fields  (A  » b).  From  (1 ) and  (3) 
and  using  the  fact  that  at  room  tempeiature  A'j  2jA'i,  wc 
obtain 

_ _J ^1  1 

ITia  + 1(^3  ' ‘*^31  )1  [I  3j  + t(cO|  - W32 )] 


Expression  (5)  applies  to  stationary  molecules  w'ith  energy 
levels  at  Ei , , and  Ey . In  a gas  sample  we  need  to  account 

for  the  Doppler  shift  of  the  transition  energies  of  individual 
molecules.  This  is  done  by  averaging  the  nonlinear  coefficient 
daffy  over  the  Maxwellian  velocity  distribution  function  with 
the  result,  for  operation  at  line  center,  that 

;u),=W,-u),  _ ~'^l  (f'7)2l(f'ii)l3(f^o).33  - 

•V^/2  — ^ (6) 

OCO31  ol 


F{x)  = e^  erfc(x),W3i  =»  wjj  »toj,,r,3  =«  Fj,  s P, 

o = \/WjM 

is  the  rms  molecular  velocity,  and  x = cri(\/2cu}3,)  is  the 
ratio  of  the  homogeneous  (spontaneous  plus  pressure)  line- 
width  r to  the  Doppler  linewidth  \/Tocj3,  /c. 

Althougli  a numerical  estimate  of  the  nonlinear  mixing  coef- 
ficient based  on  (6)  is  possible,  a safer  procedure,  and  one  that 
serves  as  a check  on  the  matrix  elements  needed  to  evaluate 
dxxT  (the  largest  coefficient  in  NHjD)  is  to  relates  it  to  the 
linear  absorption  coefficient  of  x polarized  field  at  013  = £031. 
The  latter  can  be  shown  to  be  given  (in  esu  units)  by 

731  = \At  xe^'  erfc  (x)  (7) 

where  7//  is  the  value  of  731  at  high  pressures  (cP » ocj^i) 
and  IS  given  by 


Combining  (6)  and  (7)  leads  after  some  mathematical 
manipulation  to 


PRESSURE,  tor 


« to , - u) , _ ^(7*2)1;  f 

■ gn/itoj,  U.s/ 
[2X’ 

r 2 0W31  I 


2x^F(x)- 


The  various  constants  in  (8)  are  evaluated  as  follows:  the 
matrix  element  (U;)i2  is  a function  of  the  admixture  and 
according  to  the  wave  functions  (2)  is  given  by 

<4o4o|/r3|4,4i>.  (9) 

We  obtain  the  matrix  element  (irlp.ls)  from  comparing  the 
splitting  Ej  - E\ , as  given  by  (2)  to  the  experimental  tuning 


X'Cr/ 

lig.  2.  Thcoictical  dependence  of  Ntijl)  nonlinear  coefficient  on 
picssurc  wlicn  tlic  applied  fields  are  exactly  resonant  witli  the  Stark- 
tuned  energy  levels. 


curve  of  E3  - El  versus  E^^  [14].  Tliis  yields  <i7|pj|j>  = 
1.14  X lO"'®  esu.  At  resonance  = 3570  V/cm  and  6/A  = 
0.174.  These  data  are  used  in  (9)  and  result  in 

=0.174  (i7|/j|s>  = 0.198  X 10*'®  esu. 

Tltc  saturated  absorption  7/;  and  pressure  broadening  coef- 
ficient arc  obtained  from  the  data  in  [13]  as 

7//  = 0.028  cm*' 

P/P  =271(20.1  MHz/torr). 

With  these  data  we  obtain 

= 2.31  X 10*''f7(x)esu  (10) 

C(x)=2x^^-xe*’ erfc(x)j.  (11) 

The  theoretical  dependence  of  on  pressure  (see  (10)]  is 
plotted  in  Fig.  2.  The  peak  occurs  at  P = 2 torr  and  has  a value 
of 

= 6.4  X 10*®  esu  = 2.4  X 10**’  MKS. 

A comparison  of  this  predicted  behavior  wi'h  experiment  is 
given  in  [10] . 

The  coefficient  d estimated  above  refers  to  the  generation 
of  sideband  radiation  at  coj  by  mixing  a COj  P(20)  line  with 
a mictowave  field  eoj  (at  4.1  GHz).  It  is  thus  appropriate  to 
compare  it  to  the  elcctrooptic  coefficient  u,  of  GaAs  which 
can  be  used,  alternatively,  to  generate  the  sideband  by  con- 
ventional elcctrooptic  modulation. 

Using  the  correspondence  (16,  ch.  16] 


we  have 

(ri®r)wi|  n 

We  thus  reach  tlic  conclusion  that  for  sideband  generation,  dc- 
biased  NlljD  at  /’ = 2 torr  is  comparable  to  GaAs  (which  is 
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one  of  the  best  infiared  inodiiintinn  materials).  We  must 
recopni/e,  liowever,  that  tins  large  coclt'icient  was  obtained 
by  e.sploiting  the  resonant  nature  of  the  effect.  The  penalty 
we  pay  for  this  is  that  of  reduced  bandwidth. 

III.  Conclusions 

In  conclusion,  we  have  shown  in  detail  how  Stark  admixing 
can  give  use  to  second-order  optical  nonlinearities  in  gases. 
We  have  deiived  an  expiession  for  the  coefficient  describing 
the  mixing  of  an  infrared  and  a microwave  field  in  NH2I). 
.\\ailable  absorption  data  were  used  to  obtain  a luimeiical 
estimate  for  the  mixing  and  to  describe  its  parametric  de- 
pendence. An  experimental  demonstration  of  litis  effect  is 
described  in  the  following  paper  [10) . 
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ABSTRACT 

The  paper  develops  a consistent  formalism  for  describing 
nonlir.oar  optical  mixing  and  multi  photon  processes  of  any  arbitrary 
order.  The  theory  uses  the  time-evolution  operators  of  quantum 
mechanics  and  the, related,  Feynman  diagrams. 
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I.  General  Background 

The  purpose  of  this  note  is  to  demonstrate  the  use  of  the  quantum 
mechanical  evolution  operator  formalism  and  the  closely  related  technique 

of  Feynman  diagrams  in  deriving  the  nonlinear  constants  characterizing 
different  multi  photon  nonlinear  optical  processes. 

As  an  illustration  v/e  apply  the  method  to  sum  and  frequency  genera- 
tion, two-photon  absorption  and  stimulated  Raman  scattering.  We  start  with 
a short  review  of  the  evolution  operator  formalism. ^ 

The  eigenfunction  ip(i)  of  an  atom  subjected  to  an  electromagnetic 
field  can  be  obtained  formally  by  solving  the  Schrbdinger  equation 


Ht|;  = ill  II  (1) 

or,  equivalently,  by  operating  on  the  eigenfunction  at  time  t with  the 

d 

evolution  operator  u(t,  ,t  ) according  to 

0 a 

'l^(tb)  = ip(t^)  (2) 

where  '^(tb.t^)  satisfies 

u(t.  ,t  ) 

If  H does  not  depend  on  time,  then  it  follows  from  (3)  that 


-i  Jr(t. -t  ) 
u(tb.t^)  = e " 


= ^ |m>  <m|  e b a 


(4) 


I 


m 


3 


where  r and  lm>  is  the  eigenfunction  of  H with  energy  E^,  i.e., 

H|ni>  = . 

In  the  cases  of  interest  to  us  we  take  the  Hainiltonian  as 

H(t)  = H^+v(t)  (b) 

where  is  time  independent  and  v(t)  represents  the  time  dependent  in- 
teraction of  the  atom  with  the  optical  fields. 

The  solution  of  u(t^,t^)  when  the  Hamiltonian  is  given  by  (b)  can 
be  obtained  by  a perturbation  expansion  in  powers  of  v(t).  The  result  is^^^ 


u(t^j.t^)  = u^°^(t^^,tj+u(^^(t^^,tj  + u^^^V.tJ  + -.-ul"^t^,,tJ+  ... 


L-i  % W> 

b*  a 


^ r s V 

b’  a 


b’  a^ 


^n  ^''b’‘'a^ 

(6) 


where 


u^°^(ti).ta)  = e 


IT 

u^*^^(t,  ,t)  v(t)  u^^'(t,t^)  dt 


.(2)  . , i^2 


V'^i 


(-fr)^  J J dt^dt^u^^^tj^.t^)  v(t^)  U^''^(t^.t2)  V(t2)u^''^(t2.t^) 

(7) 


.(0), 


,(0), 


a a 


u<^>  = (-  i)3 


^a  ^a 


dt^dt2dt2U^*^^{tj^,t^ ) v(t^)u^®^(t^,t2)  v(t2) 
X u^^^t2.t2)  v(t2)  u^^^t^.tg) 


where 


> t,  > > tj  > tj  ...  > 


The  manner  in  which  we  are  going  to  apply  this  formalism  may  be  per- 
haps best  illustrated  by  an  example: 
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Consider  the  Kaiiidn  scattering  process  in  which  an  atom  (or  molecule), 
subjected  to  an  optical  field  containing  frequencies  u-j  and  and  initi- 
ally in  the  ground  state  n,  absorbs  a photon  at  , while  simultaneously 
emitting  a photon  at  0)2  < while  making  a transition  to  state  k.  This 
process  is  illustrated  in  Fig.  1.  Each  scattering  event  results  in  the  ab- 
sorption of  a photon  at  and  the  emission  of  a photon  at  We  may  wish 

to  calculate  the  rate  at  which  photons  are  generated  per  unit  volume. 

This  rate  divided  by  the  ptioton  flux  (photons-m"  -sec  ) gives  us  the 
exponential  gain  constant  of  the  field  (the  Raman  gain). 

We  calculate  the  rate  at  which  an  atom  makes  the  Raman  transition 
n ->■  k by  taking  the  time  derivative  of  the  probability  that  an  atom  ini- 
tially in  the  ground  state  n will  be  found  at  time  t in  state  k.  This  prob- 
bility  is 

P^(t)  = |<ki’^(t)>i^ 

= l<klu(t,0) [n>l^  (8) 

Since  the  basic  scattering  involves  two  photons,  one  at  cj-i  and  a second  at 

(2) 

0)2,  v.'e  need  to  consider  only  the  second  order  term,  i.e.,  u^  '(t,0)  so  that 

P^,(t)  = |<k|u^^^t.O)ln>l^  (9) 

An  alternative  method  would  be  to  solve  for  the  third  order  induced 
dipole  moment  of  an  atom  subjected  to  fields  E.jexp(iWit)  and  E2exp(iW2t). 
Specifically  we  are  looking  for  an  induced  dipole  moment  which  is  propor- 
tional to 

io),  t itot  ico^t  « i(jj«t 

E^e  ' (E^e  )*  ^ = jE^p  £26  (10) 

Since  the  expectation  dipole  moment  is  given  by 


\ 
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*1 


<M.>  = <i,0(t)  I^J^.  [i,'  (t)> 


(11) 


The  desired  dipole  moment  will  result  from  terms  such  as  <1;;^  ^ ^ | 1 

I p^  1 whore 


i//"^(t)  H u^'^^t.O)  ;(o) 


since  they  involve  the  third  power  of  the  perturbation,  v(t). 


(12) 


Using  the  first  approach  the  power  per  unit  volume  generated  at 


is 


P(w2)  = 


where 


Vk  = k hit’ 

is  the  transition  rate  per  atom  from  |n>  to  |k>,  and  N is  the  density  of 
atoms. 

Using  the  second  approach,  i.e.,  the  one  leading  to  (11),  we  calculate 
the  power  P(a;2)  as 


P(w2)  = 


_ F 
2 ‘^21 


(14) 


v/here  the  bar  denotes  time  averaging. 

Either  result,  i.e.,  (13)  or  (14),  can  be  obtained  in  a quick  and 
straightforward  manner  using  Feynman  diagrams.  This  makes  it  possible  to 
limit  the  considerations  to  resonant  terms  which  dominate  the  process, 
when  such  terms  exist,  while  ignoring  the  remaining  terms. 

By  insisting  that  the  results  of  (13)  and  (14)  agree  with  each  other 
we  will  learn  how  to  include  finite  linewidths  in  the  analysis. 


G 


To  calculate  the  expectation  value  of  physical  observables  as  in  (11) 
we  need  to  obtain  an  expression  -^or  the  wavefunction  to  any  desired  order 
of  perlurbaticn.  The  optical  field  at  the  atom  site  is  taken  as 


i U3,  t 


iwot 


1 I UJ-i  C <1  I tUn  U 

E(t)  = ^ e + + c.c. 

and  the  interaction  Hamiltonian  as 


v(t)  = - y • ij-  e 


t.1  ico-j’w  ^2  ^ 


+ 2”  G 


+ c.c.) 


(15) 


(16) 


The  perturbation  is  assumed  turned  on  at  tp  = at  which  time  the  atom  is 
in  its  ground  state  n.  The  eigenfunction  at  a subsequent  time  t is 


Mt)  = ii.'^°^(t)  + t^^t:  + •••  4^("^t)  + 


(17) 


where  from  (2)  and  (6) 


\i;^'^^(t)  = u^"^(t,t^)|n> 


(18) 


Using  (7)  we  write 

,.(0), 


H 


(t-t  J 


-iu  (t-t„) 


^(^^(t)  =e  t ° |n>=e  ° |n> 

* ill 


(19) 


V^^^^t)  = -4-  f e"  v(t^)  e 


iij  ^ 

t. 


- -T^  (ti-t  ) 
'1  0 


dt-j  I n> 


which  using  (4)  becomes 

^t)  = - ^ n |m>  <m|e 


v(t^)e 


-'“n'h’V 


|n>  dt,  (20) 


Since  v(t^)  is,  according  to  (16),  a sum  of  four  terms,  the  integrand  in 
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(20)  is  made  up  of  four  terms.  One  such  term,  for  example,  resulting  from 
the  part  of  v(t)  involving  E*  oxp(-iw^t)  will  yield 

t 


m _ 


1 (w -0), ) t,  -i(jj  , lunt 

* „ mn  1 1 „ int  „ o . . i . 

^ G e e dt^  m> 


t ►-OO 
0 


2F 


i ( -w,  -'Jj  ) t 
1 n 


rmn  1 iY 


--  m> 


(21) 


where  p-j  is  the  projection  of  u along  E^.  In  (21)  we  added  the  usual  con- 
vergence factor  Y by  letting  -*  The  factor  y is  then  allowed 

to  approach  zero.  This  causes  the  integral  at  t=-«  to  be  zero.  We  will 
show  later  that  if  v/e  leave  y as  in  (21)  it  will  account  correctly  for  the 
finite  lifetime  or  linewidths  of  the  transition.  The  constant  factor 
exp(iwty)  has  been  le-^t  out  since  it  cancels  out  (through  multiplication  by 
its  complex  conjugate)  in  the  calculation  of  physical  observables.  As  men- 
tioned above  \|)^^^(t)  has  three  more  terms  each  of  a form  similar  to  (21)  but 
involving  and  These  arise  from  using  the  remaining  three  terms 

of  the  Hamiltonian  (16)  in  (20).  We  thus  have 


m 


(Pi ) 


E*  e 


i(-Wl-mn)t 


Jm> 


1 ^mn 


co„„-  0,',-  1Y 
mn  1 ' 


i(u,-co  )t 
e ' " lm> 

•^'mn^  <"1- 


+ (^2) 


£*  e 


i(-m2-Wn)t 


m> 


mn  w-  w,-  iY 
mn  c 


+ (U2) 


E2  e 


i(m2-mn)t 


m> 


mn  a)„  +cao"  lY 
mn  2 ' 


(22) 


(2) 

Using  (7)  v/e  write  the  second  order  wavefunction  ij;'  ^(t)  as 


t t, 

(-  f 

0 0 


dtidt2  e 


-iH  (t-tj  -ill  (t  -t„) 

° ^ v(tj  e ^ v(tj  e ° 2 0 


n> 


In? 


= ^--wr  I I 


) -iu)  (t-t 


m s 


; 

0 0 


" ''s>  <s|  v(t^)  e.  "'121 


tn> 


X <m|  v(t2)e 


. -*“n<V'o> 


n>  dt2dt^  (23) 


Since  v(t)  contains  four  frequency  terms  and  appears  tv/i  co  in  (23)  the 


full  integration  in  (23)  will  yield  16  terms.  A typical  term,  as  an  ex- 

^2  and  uu 


ample  involving  the  use  of  -^1^*  e at  t,  and  -y^E^  e ^ at  t,  is 


12  E*  E*  (vJ-i)mn(^'2^sm  ^ 


i (-to^-o.'|-a32)t 


j]  t <^2 


m s 


S> 


(24) 


A convenient  way  to  represent  (24)  is  through  the  use  of  a Feynman 

diagram  as  shown  in  Fig.  2.  Time  increases  from  the  bottom  to  the  top. 

Each  solid  line  segment  represents  an  eigenstate.  The  atom  starts  at  t^  in 

state  n,  scatters  at  t2  into  state  ni  by  absorbing  a photon  at  w-j . This 

— 1 to*,  t 

scattering  is  accounted  for  by  the  factor  (y^)|^^^E^  e 
(24).  The  next  "scattering"  is  at  t^  and  involves  an  absorption  of  a 
photon  at  ^2-  A negative  frequency  denotes  absorption  and  is  represented 
by  an  arrow  terminating  at  a corner,  while  an  arrow  starting  at  a corner 
denotes  the  emission  of  a photon.  A case  where  the  transition  from  state 
n to  s involves  the  absorption  of  a photon  at  and  the  emission  of  a 
photon  at  ci)2.  as  an  example,  is  shown  in  Fig.  2b. 
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The  corresponding  contribution’  to  (t)  can  be  written  by  inspection 

l(-tj  -a),+(D„)t 

I (2)  = V V ^ , 


We  note  that  each  scattering,  i.e.,  each  corner  in  the  diagram  contributes 

one  factor  to  the  denominator  and  the  factor  is  equal  to  the  enemy  in 

units  of  fi  of  the  atom  and  field  after  the  scatterino  minus  the  initial 

(t  = t^)  energy.  The  second  factor  in  the  denominator  in  (25),  as  an 

example,  is  obtained  from  Wp+tj2-  (w-i+w^).  The  remaining  fourteen  wavefunc- 

(?) 

tions  are  obtained  by  taking  all  tlie  possible  permutations  of  o)  _ , 

(2)  -w. 

i»j=l>2  where  our  corivention  is  such  that  _ . for  example,  corre- 

^i 

sponds  to  emitting  an  photon  at  t^  and  absorbing  an  Wj  photon  at  t-j 

(t,  >to).  We  note  that  is  not  equal  to 

'12  ^ 

Using  the  diagram  technique  we  can  write  the  wavefunction  to  any 

order  of  perturbation.  As  an  example  consider  the  process  in  which  an 

atom  makes  a transition  from  state  n to  state  s v;hile  absorbing  one  photon 

at  tu-j,  one  photon  at  ^2  and  emitting  a photon  at  ai^.  One  diagram  describ- 

(3) 

ing  this  process  is  ’A  ' as  shown  in  Fig.  2c.  We  obtain  by  inspec- 

^ 2 ^ ^2  ^ 3 


” (0  j “ (1)2  • 2 ” 


\\\  ^3  ^^^1  ^mn^^2^km^^^3^sk  ^ 

m k s 8h^  ((0  ) (m.  -Wo ) ( 


i(-to^-a).|  -02+^2)^ 


(3)  3 

The  total  number  of  \lr  ^ , combinations  is  6 = 216.  The  tremendous 

-W.  ,-lOj 

advantage  of  the  diagrammatic  representation  is  that  it  enables  us  to  single 


1 


out  and  then  treat  very  simply  the  terms  which  dominate  in  any  given  physical 
situation.  This  will  become  clearer  in  the  examples  whicli  follow. 
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II.  Froquency  Addi IJ on 

In  this  section  v/e  will  calculate  the  nonlinear  susceptibility 


'^ijk 


which  describes  via  the  relation 


(j3-co^+cj2 

= ^ijk 


(27) 


The  induced  polarization  along  the  i direction 


/ au'-Wi'ioj..  iw^t 

P^Ct)  = Re|p/  ^ e ^ 


ICO-j  t 

due  to  input  fields  Re  E.jj  e jand  Re 


\hk  ^ 


i a'2t 


(28) 


polarized  along  j and 


k,  respectively. 

The  total  atomic  dipole  moment  is  obtained  from 


<y^-(t)>  = <ijj(t)  lii.  |4y(t)> 


(29) 


so  that  contributions  containing  the  electric  field  to  second  power  as 
in  (27)  must  come  from 

<yj^^(t)>  = iy.j  |y.j  ^ |u.j  (30) 


Each  term  in  (30)  contributes  two  terms  to  the  desired  polarization--one 
in  which  E^  scatters  first  and  then  E^*  and  vice  versa. 

Consider  as  an  example  the  terms  due  to  ^ 1 y.j  1 4'^^  : 

i(aj2+u)n)t  i(a).|-(a^)t 

Im> 


i (aj2+u)^)t 

^ ^ ^^'''2^^'k^sn  ^ -0)2+  iy  1^‘il^l^^j^ 

s m 4ti  sn  2 ' 


n,n  ly 


' I I 


s m 4h 


^1^2  ^i(a)^^ta2)t 
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Note  that  any  constant  prcMii/lti plying  a v/avefunction  to  the  left  of 
is  replaced  by  its  complex  conjugate.  Collecting  all  the  terms  con- 
tributing to  yields 

y y [ ^^k^sn^^i  ^sm^^j  ^mn  ^ ^^j^sn^'S'  |sm^^k^mn 


+ ^^j^nir/^VJsni^^i  ^is ^ ''  k'mn  ’'"i  'sni''~1  'ns 

T^a^^-iYKco^n+oj^+oj^-iY)  TI>^n'^‘^32-i7Tro’s,+a;^+to2-'"^' 


“sn  “'I  “^2 


^ ^ns  ^ns  1 

A more  conventional  derivation  is  given  in  References  (2),  (3)  and  (4). 
III.  Two-Photon  Absorption 

Here  we  will  apply  our  formalism  to  derive  ti'.e  transition  rote  for 
absorption  of  two  photons--one  at  to^  and  a second  at  aj2--by  an  atom.  By 
allowing  we  will  get  the  familiar  expression  for  two-photon  absorp- 

tion coefficient. 

We  start  by  writing  the  second  order  wavefunction  i/)^^^(t)  correspond- 
ing to  an  atom  which  at  t=0  is  at  state  n and  which  interacts  with  the 

radiation  fieldi  Eq.  (15),  consisting  of  fields  at  and  We  will 

(2)  (2) 

assume  that  the  largest  contribution  to  \J;  comes  from  ' corres- 

"co-j  I " ^2 

ponding  to  the  diagram  of  Fig.  3.  We  write  by  inspection 

-iw^t  r i(io^^-w^-a)2)t  , 

|s>  (31) 


(t)  = y y (^i^)  ^ 


m s 


The  -1  term  is  due  to  the  fact  that  here  we  integrated  from  0 to  t.  Let 

us  assume  that  for  some  level  m,  w = oj,  so  the  term  m dominates  and  we 

mn  1 


I rewrite  (31)  as 

t 

I 
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r -iwt  i(w  -to,-Wp)t 

v,(2)(t) . ‘1^  I — !_ ; 

5 I Stl  1 1 


The  probability  of  finding  the  atom  in  some  state  k at  time  t is 

p 

|<k|iij(t)>|  so  that  the  transition  probability  due  to  tv/o-photon  absorp- 


tion is 


|<kk^^^t)>|2 

16fi^  (%in-“l)^ 


We  note  that  transitions  occur  to  the  state  k which  conserves  energy, 
i.e.,  to  that  state  where  = 'h((0i+i02) . If  the  normalized  lineshape 

function  for  the  transition  n ->•  k is  the  average  value  of 

P|,  is  obtained  by  multiplying  (34)  by  g(tJJ[,p,)  arid  integrating  from  -«•  to 
00^  Using 

sin^  xt 

lim  ^ — 5-  2-iTt  6(x)  (35) 

t (x/2)‘^ 

we  obtain  for  the  transition  rate 


W = — P 
”n  ->k  dt  '^k 


_ g(^kn^^l'^^2^  (3^) 

The  special  case  of  one  frequency  two-photon  absorption  where  cji= 

(5) 

has  been  considered  in  some  detail  . We  may  obtain  the  absorption  coef- 
ficient for  this  case  by  equating  the  change  in  the  intensity  I to 
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the  number  of  transitions  per  unit  time  per  unit  volume. 


dl 


to 


“cfT 


2)ito  W 


L (N, -N  ) 
n ► k i;  n 


(37) 


2 

where  is  tlie  intensity  (watts/ni  ) at  to  and  N 
tion  densities  of  level  k and  n,  respectively, 
obtai n 


and  are  the  popula- 
Usinq  1^^  =^|E,|^/2  we 


“two 

photon 


mo(N^-N'k)(ni)^,„(l.p^,/i^  p(coi-,n-2to)l^ 

-h^(to  -co)^ 

' mn 


(38) 


IV.  Raman  Processes 

As  the  most  detailed  example  of  the  diagram  formalism  we  take  the 
case  of  Raman  processes  in  which  an  atom  initially  in  the  ground  state  n 
absorbs  a photon  at  to-j,  emits  one  at  while  making  a transition  to  the 
state  k.  The  situation  is  depicted  in  Fig.  4.  The  situation  is  identi- 
cal to  that  considered  in  the  previous  section,  except  that  here  the 
second  photon  at  [O2  is  emitted  rather  than  absorbed.  This  causes  the  0^2 
arrow  in  the  Feynman  diagram  of  Fig.  4 to  have  the  opposite  sense  to  that 
in  Fig.  3.  We  can  use  directly  the  result  of  Section  III  and  merely  re- 
verse the  sign  of  0:2-  Using  (34)  we  obtain  for  the  transition  probability 


Pk(t)  - — -T- 2TTt  6((ak„-u^+a>2) 

let,  (m^^-a)^) 


(39) 


In  (39)  it  was  assumed  that  one  term  m dominates,  so  that  the  sunriation 
over  m is  omitted.  As  in  (38)  we  introduce  the  normalized  Raman  line- 
shape  function  g(mkj^)  such  that  g(tL'kj^)dUk^  is  the  probability  of  having 
within  the  interval  dcj^^  . Multiplying  (39)  by  g(c)^„)dol^„  and  inte- 
grating leads  to 


kn'^'kn 
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which,  using  (42)  gives 

9r1  ”3'  2 ' ; a)^+(D2) 


Next  we  will  describe  the  same  Raman  scattering  leading  to  (43)  by 
a different  approach.  Instead  of  using  the  concept  of  a transition  rate 
we  will  calculate  the  induced  dipole  moment  at  of  an  atom  interacting 
with  the  and  fields.  To  make  sure  that  v/e  will  treat  the  problem 
to  the  same  order  of  interaction  we  look  for  a dipole  moment  cubic  in  the 
applied  field  amplitudes. 
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<M.>  - E,E*E2  e 


i((i)l-u)i+  (1)2  )t- 


= lEJ' 


E2e 


1tL)2t 


Such  a contribution  to  which  involves  the  third  power  of  E can  come 
from  combinations  such  as  ^ |u  I U J An  examination 

of  (22)  and  (25)  shows  that  the  largest  dipole  moment  will  arise  from  the 
use  of 


(44) 


<p^.>  = ^ (t)  I p + c.c. 

where  taking  the  direction  of  E^  as  j and  as  £ 

- i ‘zr  (10  -n  -0.,)  ° 1“’ 

m ' mn  'm  r 

i (-a),+u„-(jj  )t 

g s 4h"  (tx'g^-a)^-iY)(aJ^^-w^+W2-iY/ 

so  that  using  (44), 


<I'1>  ” -V  I n EiE^Ej  J“2\ 

W mgs  (“mn-“l'<“gn-“l><“sn-“r“2'*^> 


s> 


(45) 


(46) 


v/here  we  dropped  the  term  iY  in  the  first  two  factors  in  the  denominators. 

These  must  be  reintroduced  if  operation  is  such  that  ui  - to,  or  oj  ='  to,. 

mn  1 gn  1 

We  will  assume  that  the  summation  over  m is  dominated  by  one  term 


s = k 


‘"mn  " ‘"r 

Also, 

‘"kn  = ‘"r 

• (O2.  1 

^ 1 

<p^>  = 

’ eh-^ 

(to  -to 
' mn 

ii02t 


EiE*  E2-  c.c. 


(48) 


We  define  a nonlinear  Raman  susceptibility  by 
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P,  ^ = N<„,>  ■ Ej 


(49) 

From  (48) 


(cOq  ) 

where  P is  the  complex  amplitude  of  the  polari7ation  at 
and  (49)  v/e  obtain 


~ 3 2 (>>0) 

The  average  power  per  unit  volume  generated  at  is 

(0  G f\  r\ 

P = " 1^2!  Xr  (51) 

where  Xr  = Xr  - iXp-  (52) 

Ar.  examination  of  (50)  shows  that  when  y = 0 Xr  Is  real  and  p - 0.  The 
convergei^ce  factor  y 'Ti  the  Feynman  formalism  is  thus  fundamentally  related 
to  the  problem  of  power  dissipation.  Let  us  for  a moment  consider  y as  a 
parameter  and  calculate  the  power  p.  We  obtain 


I ^ I < I rm  I i iiir*.  i\il  I 


(2,  .2- 


and  from  (51)  (putting  i=2  since  only  the  polarization  component  along 
contributes  to  absorption)  ^ 2 

. _ ‘^'2  ic  |2ir  |2  ^'^l^nin  ^^2^mk  ^ 

P~'~::3i^l'  '^2  ? 2 7 — 


N 

wL' 

^‘^mn 

-N 

.^3 

4h 

(o)  <-0) 

' mn 

This  power  must  be  the  same  as  that  calculated  using  the  transition  rate 
approach  leading  to  (41).  The  latter  leads  to 
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p = NW„  ^.a>2 
Nnuj  |E,|  lEjl 


8h' 


(i-5) 


Equating  (5^)  to  (5D)  and  taking  y^=U2 


Y/tt 


^‘"'kn"^^'r“2^^  + Y 


'2  2^  " 9[^[,p“  (<^l"W2)l 


(9G) 


We  thus  obtained  an  explicit  expression  for  the  transition  lineshape  func- 
tion g[c0|^^-(a)-j-t)j2)].  We  now  find  that  the  result  of  the  coherent  treatment 

in  which  the  problem  is  treated  by  the  induced  dipole  agrees  exactly 
with  the  transition  rate  approach  if  we  merely  regard  y as  the  Lorentzian 
width  of  the  k->n  transition. 


V.  Two-Photon  Raman  Processes 

As  the  last  application  of  the  evolution  operator  techniques  we 
consider  the  problem  of  two-photon  Raman  transitions  of  the  type  consid- 
ered in  Fig.  5.  These  processes  were  described  by  Yatsiv  et  al.^^^.  We 
shall  expand  this  discussion  and  show  how  these  processes  can  be  de- 
scribed in  the  formalism  of  nonlinear  optics. 

In  the  first  of  these  processes  {5a)  one  photon  at  co^  and  one 
photon  at  0)2  are  absorbed  and  a photon  at  is  emitted,  while  the  atom 
(molecule)  makes  a transition  from  the  initial  state  n to  some  final 
state  k.  In  the  second  process  one  photon  (to^)  is  absorbed  while  photons 
0)2  and  o)^  are  emitted  during  the  same  transition. 
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The  analysis  of  both  these  processes  is  similar  and  the  results 
for  5b  can  be  obtained  from  those  of  5a  by  reversing  the  sifjn  of  ^2- 

Using  the  formalism  of  nonlinear  optics,  the  scattering  events  of 


Fig.  5 can  be  represented  by 


u)3=m^-<a^+U2 -02+0,3 


'^o^R2 


(Oo=0, -0,+(jJo-Oo+Oo)  O o 

J I I C C S I r \ \ r \ ^ 


so  that  we  may  expect  an  amplification  (or  attenuation)  at  03  with  an 

2 2 

exponential  constant  proportional  to  |E.||  IE2I  in  analogy  with  the  one- 
photon  Raman  process  (43).  Alternatively,  we  can,  as  in  Section  IV, 
obtain  the  gain  by  calculating  the  transition  rate  per  atom  for  the 
process  illustrated  in  Fig.  5a. 

By  direct  analogy  with  the  calculation  leading  to  (39)  we  obtain  for 
the  probability  P|^(t)  of  finding  an  atom  in  state  k 

(t)>|  = :T,  ;2  2Tit6(a)|^^-w^-t^2'^‘'^3' 

<"sn-“r“2>  (jj. 


SO  that  the  transition  rate  ^ is 

ttI  ^ I21 3(01  )^n(P2)ms^P3^sk"l"2"3^  1^2 

^'n->^k  3 3,  32,  32 

4h  c r.  (co^^-m^) 


g(a,kn=^l+(^2’‘^3^ 


The  power  per  unit  volume  generated  at  0,3  is  thus 


P(w3)  - '(,0,3  - 91^2^3 

where  the  two-photon  Raman  gain  g|^2 

V"'3(Pl^m(P2\L(P3’sk  h^2  , 

^R2  “ TTTTX  72~,  ^(‘*’kn'‘"’l  ‘‘'2’“3^ 
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and  is  thus  proportional  to  the  product  of  intensities  at  and 

We  note  by  comparison  to  the  expression  (43)  for  the  one-photon  (i.e., 
conventional)  Raman  gain  that 


9r2  ""  ^R1  ^2 


2h  ceCto^^-w-j-w^) 


(6?) 


in 

If  we  assume  that  I2  10  watt/cm,  '^10  , i.e.,  near  resonant 

-10 

tuning,  y n.  10  MKS,  v/e  obtain 

^R2  ~ ^R1 

so  that,  using  intense  laser  sources  tuned  near  atomic  resonances,  it  may  bo 
possible  to  obtain  two-photon  Raman  gain  approaching  those  of  normal  Raman 
processes . 

Unlike  frequency  addition  and  multiplication,  this  process  does  not 
require  phase  matcliing.  This  is  due  to  the  fact  that  if  the  spatial  phase 
factors  are  included  in  (57)  they  will  cancel  out,  since  each  field  ampli- 
tude is  multiplied  by  its  complex  conjugate  except  for  E,.  The  induced 
polarization  P thus  has  the  same  spatial  phase  as  and  the  interac- 
tion is  spatially  cumulative.  If  the  interaction  region  is  provided  with 


feedback,  osci  1 lation  at  niay  result. 

To  describe  the  same  process  using  the  concept  of  nonlinear  Raman 
susceptibility,  we  proceed  to  derive  the  dipole  moment  per  atom  <y^>  at 
induced  in  the  presence  of  fields  at  w-j , W2.  and  Resonant  contributions 

to  the  fifth  order  dipole  moment  arise  from 

<y.>  |y. 

^-a)i,-u)2'  1'  -Ui,-t02,t 
From  (25)  and  (26)  we  obtain 

r*  r*  « 1 2 n , 


> + c.c. 


1 .2  Et  e 


-a)^-(j02-iY) 
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s> 


(63) 


1 1 1 
m k 


i (-io,-aj^+w^-(j  )t 

<“lln‘“r’’'><"rml'"r“2-'''>‘"krr“r“2*“3-*''> 


SO  that 


<u^>  - ( — 3-)  I 

32h  gsi-nik 


|E-|I  Ugl  ^3^^‘i^an^^‘2^sq^'^'!^nil^'‘?^£m^^3\nk^^‘i^sp'  ^ 

(u)^^ -c^l  - i Y ) ( -''^r‘^'2 - ' ^ ^""kn ""^r‘"'2^  V ’ ^ -"^1  ^ ^ ‘"sn V ' ^ 


- — + c.  c. 


Let  us  assume  that  to^  and  co^  ^ne  adjusted  so  that  for  some  levels  P- 
and  m,  to.  = to,  while  to  = to,+(o„.  Also,  that  a level  k exists  such 
that  t0|^^  = t0i+t02-t02.  Keeping  only  the  resonant  term  of  (65)  we  obtain 
9 ? 97 

<p  > ~ 1 1^1 1 1^2!^  ^3^^1-g.n^^'2^icr'^^3^niiJ‘^kn~^‘^l'^^'2"^'3'^'*>-n^^  , ^ 

(66) 

If  we  relate  the  complex  amplitude  P^(t03)  of  the  induced  polarization 
N<vi^>  to  the  complex  field  amplitudes  via 

^i  ^ ^o^R2l^ll  I ^2 1 ^3  ’ ^R2  ""  ^R2  “ ’^R2 


we  obtain 


N ^^1  ^£.n^^2^^m^^3^mk^''i  ^mk 


N ^^l^L^^2^L^'‘3^mk^^i  ^mk  ^pn 


(68) 

where  2Y^j  is  the  full  width  at  half  maximum  of  the  i -♦  j transition. 
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It  follows  jntnediately  that  the  presence  of  a negative  X|^2  causes 


an  amplification  at  with  an  exponential  gain  cons  tan 


t(7) 


9r2  ""  " cn^  1^2! 

N^‘^3(h)L(''2)L(^‘3)mk  ^V3^1^2 


(69) 


(70) 

2 

where  we  used  I.  = ce.lEj  /2n..  The  result  (70)  for  Od,  is  identical 

1 1 1 1 * Kc 

with  (61)  once  we  associate  the  second  factor  in  (70)  with  the  normalized 
lineshape  function  g[a)|^^-(tii)^ +012-0)2) ]. 

We  thus  find  that  the  gain  exercised  by  the  wave  proportional 

to  the  product  I1I2  of  the  intensities  at  oi^  and  012.  so  that  a stimulated 
Raman  emiosion  at  oi^  may  be  expected  at  some  critical  value  of  Iil2- 

The  analysis  of  the  two  photon  (emissive)  Raman  process  shown  in 
Fig.  5b  is  similar.  All  we  need  to  do  to  obtain  the  gain  at  oi^  is  to 
replace  0)2  by  -oi2*  A fundamental  difference  between  the  two  processes, 
however,  would  be  revealed  had  v/e  considered  the  consequence  of  the  quan- 
tized nature  of  the  field  operators.  The  temporal  rate  of  the  photon  number 
at  oi^  corresponding  to  process  5a  is 


^=6^0^02(03+1)  (71) 

while  that  of  5b  is 

dn.,  dop 

Hr  " dF“  = C2n^(n2+l)(n3+l) 

where  n^  is  the  number  of  photons  at  to^ , and  C2  are  rate  constants. 

It  follows  immediately  that  the  stimulated  emission  of  photons  in  5a 


\ 


requires  the  s>.’Tiultaneous  presence  of  the  fields  at  and  while  in 
5b  stimulation  can  take  place  with  a single  input  at  , since  dn^/dz  ' 0 
when  n^  = n^  = 0. 

We  may  also  note  that  the  emissive  two-photon  Raman  process  of  5b 
can  be  used  to  generate  long  wavelength  radiation  at 

To  summarize,  the  formalism  of  time  evolution  operators  and  the 
related  technique  of  Feynman  diagrams  was  applied  to  a variety  of  multi- 
photon  processes  in  nonlinear  optics.  The  formalism  affords  relatively 
direct  and  orderly  treatment  of  complicated  high  order  processes. 
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